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Abstract 

A simple burning or epidemic type of algorithm is developed in order to test whether any loops 
in percolation clusters link a fixed reference loop, a problem considered recently by Gliozzi, Lottini, 
Panero, and Rago in the context of gauge theory. We test our algorithm at criticality in both 2d, 
where the behavior agrees with a theoretical prediction, and in 3d. 
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I. INTRODUCTION 

Recently, Gliozzi et al. have studied percolation in the context of gauge theory. They 
considered the question of whether closed paths in three-dimensional percolation clusters 
are linked topologically to given closed loops, the so-called Wilson loops. Studying this 
problem for rectangular planar loops, and in comparison to percolation in three-dimensional 
slabs which they relate to the problem of deconfmement, the authors find a new universal 
amplitude ratio. This work provides an example where the percolation model possesses 
connections to fundamental problems in theoretical particle physics. 

Gliozzi et al.'s numerical results for rectangular loops of dimensions RxT confirmed the 
expected behavior for p ^ p c jjj 



(W(R,T)) = Ce-^'^R^J^j^ , (1) 

where (W(R,T)) is the average probability that there is no path in any cluster linked to 
the Wilson loop, C, P, and a are constants that depend upon the percolation probability 
p, and i] is the Dedekind function 77(f) = g 1 / 24 fj^ =1 (l — q n ) with q = e 2l7TT . When p < p c , 
one expects a = because the linking probability should depend only upon the perimeter 
of the loop, while for p > p c , the linking probability is expected to decay exponentially with 
the area of the loop (a > 0). Taking p somewhat above p c , Gliozzi et al. found that the 
dependence of a upon p behaves as 

a = S(p-p c ) 2 » , (2) 

similar to the behavior of a surface tension, where v « 0.8765 is the correlation-length 
exponent of 3-d percolation □ and S is a constant. The y a.so defined the portion 
threshold pg for slabs of thickness I in the range 3-8. One expects t ~ £,(pe) ~ T~ 1 (p£— p c )~ u 
where T c is a constant, and indeed they find 1/ l {i^Ja{p i )) ~ T c /\/S is a universal amplitude 
ratio with a value of about 1.50. 

In this note, we discuss two points related to the work of Gliozzi et al.: (1) We describe an 

n 

epidemic or burning 5] type of algorithm that may be simpler than the algorithm described 
by Gliozzi et al., and (2) we apply it to study the linking probability exactly at p c (a point 
that Gliozzi et al. did not consider) for 2-d and 3-d systems. Note that Eq. (Q) and its 2-d 
analog are not necessarily expected to be valid at p c . 



II. ALGORITHM 



Gliozzi et al. describe an algorithm that involves successive removal of dangling ends and 
reduction to an auxiliary graph that represents the connections between clusters on either 
side of the flat region enclosed by the loop. This graph is used to determine whether a 
cluster is linked to the loop. 

Here we describe a cluster burning type of algorithm that accomplishes the same test. 
As in Ref. we consider the loop 7 to be on the dual lattice, so effectively the problem is 
to find if there are clusters that simultaneously pass through the plane E of vertical bonds 
enclosed by 7 and through bonds in the same plane outside of 7. We are thinking of a simple 
cubic lattice to be specific. 

To begin the process, all of the sites are set to the "unvisited" state and bonds to the 
"undetermined" state. Then we pick one of the (unvisited) sites S directly above E, and label 
that site as "visited" with an arbitrary index n. We check the six bonds that emanate from 
S; the undetermined bonds are made "occupied" with probability p and "vacant" otherwise. 
For the bonds that are occupied, we check the adjacent site; if that site is unvisited, we label 
it as visited (with a value of the label described below) and put its coordinates on a queue 
for future checking. After finishing checking all bonds connected to the site being studied, 
we consider the next site on the queue, continuing this process until the queue is empty. 
This is the normal burning or epidemic type of algorithm to identify a cluster connected to 
a site in bond percolation; here we also decide whether a bond is occupied or not as we go 
along. We repeat this process for all remaining unvisited sites in the plane above £. 

What we now do differently for the loop-linking problem is that we assign an index n to 
each visited site in a cluster. When we transverse one of the occupied bonds that intersects 
E, we increment n by one when going downward or decrement it by one when going upwards. 
In this way, every site of the cluster will be labelled by n, n ± 1 (if a path of the growing 
cluster goes once through E), n ± 2 (if a path of the growing cluster winds twice through 
E), etc. 

Now, if during the growing process a new bond is found to connect two visited sites with 
different labels £, then the cluster must have wrapped around 7 and is therefore linked to 
it. 

An occupied bond of course will not connect to sites of two different clusters, by definition, 
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FIG. 1: Logarithm of (W(R)) (= the probability that no cluster encircles just one of the two 
points), vs. the logarithm of the points' separation R. 

so therefore one does not have to worry about interference between clusters in this algorithm. 

Here we have not dealt with the system boundaries. If open boundaries are used, care 
must be taken so that the boundary bonds are not mistaken for wrap-arounds. To eliminate 
this problem and to lessen the effects of the boundary, we considered periodic boundary 
conditions in all directions. This makes the problem slightly different, because periodic 
wraparounds through the Wilson loop will also contribute to linking events, but if the 
lattice dimension L 3> R, this difference should not be too significant. Indeed, a moment's 
reflection shows that it is very unlikely that there will be wraparound without linking, since 
a cluster that wraps around the lattice it is usually a ubiquitous one and most likely will 
link the Wilson loop also. 

The idea of adding a label to sites in percolation to test for a crossing criterion has been 
used previously in relation to wrapping a periodic system in a given direction jfllll- 



III. LINKING IN TWO DIMENSIONS AT CRITIC ALITY 

For the 2-d system, the question that is studied is whether there exists a closed path in 
a cluster that encircles one (but not both) of two points on the dual lattice, separated by a 
distance R. In this case we can make a simple theoretical prediction for (W(R)), since the 



condition of a path not encircling either of the two points individually is equivalent to the 
existence of a continuous path between the two points on the dual lattice — that is, a cluster 
that connects the two points. The density drop-off from any point on a given cluster goes as 
r D ~ d where D is the fractal dimension and d is the spatial (Euclidean) dimension. To find 
the probability that two given points are connected, the above factor must be multiplied by 
the probability the size of a cluster connected to one of the points is at least large enough to 
reach the other point. At criticality, the probability that the number of sites connected to 
a point is equal to or greater than s is given by P> s ~ s 2 ~ T where r is the size distribution 
exponent, and this implies that the probability that the radius is greater than or equal to r 
is given by P> r ~ r D ^~ T \ s j nce s ^ r D _ Then, D y the hyperscaling relation d/D = r — 1, we 
have P> r ~ r D ~ d _ Thus, the net probability that two points separated by r are connected 
by a cluster at p c is given by P{r) ~ r 2 ( D - d )^ implying that 

(W{R)) ~ iT 2(d - D) = exp(-2(d - D) In R) (3) 

In d = 2, D = 91/48 and 2(d - D) » 0.208. 

We carried out simulations for this system using the algorithm described above. We 
considered bond percolation on the square lattice at p = p c = 1/2, on a system with a 
square boundary of dimensions 1024 x 1024, and considered separations of the two points 
ranging between 10 and 100. Fig. 1 shows the results for a plot of InW vs. lni?, for a 
relatively small number of runs (100000 each). The slope is about —0.22, consistent with 
the theoretical prediction above. To make this work more precise, one would have to consider 
different size systems to study the finite size corrections, and perhaps also consider systems 
with open boundary conditions for comparison. 

IV. LINKING IN THREE DIMENSIONS AT CRITICALITY 

For the 3-d problem, we consider bond percolation on the simple cubic lattice, and take 
p = 0.2488126, which is an estimate for p c believed to be within about 5 ■ 10~ 7 of the 
v alue B. We — a latte of size 128 x 128 x 128, and square Wilson loops 
containing R x R vertical bonds, with R = 5, 11, 21, • ■ ■ ,111. Between 300,000 (smaller R) 
and 13,000,000 samples (larger R) were generated for the different values of R. 

In Fig. 2, the lower curve represents \n(W(R,R)) as a function of R. The data shows 



5 




FIG. 2: 3-d data. Upper (o): \n{R~ 1 l A {W{R, R))) vs. R. Lower (o): \n(W(R, R)) vs. R, which 
shows a good fit to a straight line for 5 < R < 111. The equation of the linear fit is given. 

quite linear behavior up to R = 111. Evidently, the perimeter term proportional to P in 
Eq. (JTJ) dominates; as expected, there is no term proportional to the area. 

We do not see evidence of the R 1 ^ term in Eq. ((TJ. The data marked by circles in Fig. 
2 represents ln(R~ 1 ^ 4: (W(R, R))) vs. R, and the fit to a straight line is much worse than for 
the case without the factor of i? -1 / 4 . This factor would show up as a logarithmic term in 
the plot of Fig. 2. 

To check further for logarithmic terms, we plot in Fig. 3 the quantity InW + 0.0765-R, 
where the constant 0.0765 was adjusted to get the best horizontal region in the center, along 
with general monotonic behavior. We see two corrections to the straight line: for small R, 
there is a small decrease, which could be fit to a very small logarithmic term, — 0.03 In R, 
much smaller than the —(1/4) In R term that would appear for p > p c according to Eq. (fljl. 
The coefficient is so small that the existence of a logarithmic term seems unlikely. 

For large R, the deviations from linearity are also small, which is surprising given that 
we went up to R = 111 in a system of size L = 128. When R approaches L it should be 
more difficult to create a linking clusters, because there is a smaller region external to the 
loop (also taking into account the periodic b.c), but this may be balanced by the increase 
in vertical wraparounds (which our algorithm takes to be linkages) through the periodic 
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FIG. 3: In IF + 0.0765.R vs. R for the 3-d data, showing deviations from simple exponential 
behavior. Error bars show two standard deviations of statistical error. 



boundary conditions. Note that at R = 111, W = 0.00028617 — that is, only 3434 of the 
12,000,000 samples did not have a linkage to the Wilson loop. 

Thus, ignoring the possible small logarithmic term, the data for 3-d (in the central range) 
yields In W = -0.07Q5R + 0.25, implying 



or P = 0.0382 for bond percolation on the simple cubic lattice at criticality. Note the linear 
fit for W above is somewhat different than that given in Fig. 2, which is just a simple linear 
fit through all the data points. 

V. CONCLUSIONS 

We see that a simple burning type of algorithm can be constructed to find the loop 
linking probability studied by Gliozzi et al. We have checked it in two dimensions at the 
critical threshold, where the linking probability is known exactly by virtue of its being dual 
to the two point probability. Of course, in 2d one can easily simulate the dual problem of 
connecting the two points. However, in 3d, where a dual-lattice procedure would be much 
more complicated, a direct determination is preferable and the algorithm presented here is 



(W(R,R)) = 1.28 e 



0.0382(2/?) 



(4) 
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efficient and simple. For 3d, we find a simple exponential relation between (W(R,R)) and 
R reflecting a perimeter effect; there is no evidence of a logarithmic correction implied by 
the i? 1//4 term in Eq. ((TJ) (which is not necessarily expected to be valid at p c ) or as suggested 
by the behavior in 2d. 
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